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Labyrinthine patterns arise in two-dimensional physical systems submitted to competing interac¬ 
tions, ranging from the fields of solid-state physics to hydrodynamics. For systems of interacting 
particles, labyrinthine and stripe phases were studied in the context of colloidal particles confined 
into a monolayer, both numerically by means of Monte Carlo simulations and experimentally using 
superparamagnetic particles. Here we report an experimental observation of a labyrinthine phase 
in an out-of-equilibrium system constituted of macroscopic particles. Once sufficiently magnetized, 
they organize into short chains of particles in contact and randomly orientated. We characterize the 
transition from a granular gas state towards a solid labyrinthine phase, as a function of the ratio 
of the interaction strength to the kinetic agitation. Spatial local structure is analyzed by means of 
an accurate particle tracking. Moreover, we explain the formation of these chains using a simple 
model. 


I. INTRODUCTION 

Labyrinthine phases are intriguing two-dimensional 
patterns occurring in various domains of physics, in equi¬ 
librium and out-of-equilibrium situations. Two distinct 
phases form at small-scale well separated stripes, which 
are themselves entangled, leading to a complex large- 
scale pattern. These shapes were experimentally ob¬ 
tained for extremely varied two-dimensional (2D) sys¬ 
tems ranging from ferrimagnetic garnet films [T] in con¬ 
densed matter and Langmuir monolayer [2] in soft mat¬ 
ter, to granular fluid suspension in which air pene¬ 
trates [1111, ferrofluid drop H and biphasic ferrofluid- 
oil layer d in fluids mechanics, and chemical reaction- 
diffusion systems [7] . The common denominator of these 
systems is the competition between long-range repul¬ 
sion and short-range attraction, which leads to the phe¬ 
nomenology of modulated phases [5]- Moreover, a wide 
range of ordering effects that lead to different patterns 
can also be related to the competition between interac¬ 
tions and geometrical frustration, as specifically shown 
for magnetic thin films d- 

By analogy with the phenomenology of these con¬ 
tinuous systems, labyrinthine and stripe phases have 
been introduced for systems of particles. In particular 
in the context of colloidal monolayers, several Monte- 
Carlo simulations dOHS] and one molecular dynamics 
simulation m have been performed. It was shown 
that, to observe stripes and labyrinthine phases, a long- 
range repulsive potential is needed, together with a 
short-range attraction which can be replaced by a core- 
softened potential da E]. Tuning geometrical frus¬ 
tration in non-interacting colloidal monolayers da US] 
leads also to stripe phases. Moreover the only exper¬ 
imental observation of a labyrinthine phase in a col¬ 
loidal system was obtained using superparamagnetic col¬ 
loids under a magnetic field, inducing dipolar interac¬ 
tions m- Labyrinthine phases were indeed found as 
equilibrium states at high enough density of micrometric 
spheres in agreement with dedicated Monte-Carlo simu¬ 
lations dSlIIS]- In contrast, similar labyrinthine or stripe 
phases have not been described in a macroscopic and out- 
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FIG. 1. Experimental setup. The 5000 chromed steel spher¬ 
ical particles (diameter a = 1 mm and mass m = 4.07 mg) 
are vertically vibrated (acceleration F = 21.9 m.s“^) inside a 
horizontal, square aluminium cell (area 5'o = 9 x 9cm^) with 
a rough bottom plate and a polycarbonate top lid (gap size 
h = 1.5a). In presence of a vertical magnetic field B, parti¬ 
cles repel each other with a force Fm- The region of interest 
is of area S = 5.7 x 5.7 cm^. 


of-equilibrium system whose particles can be individually 
identified. 

Here, we report the observation of such a labyrinthine 
phase. A monolayer of soft ferromagnetic spheres of 
1 mm in diameter is vibrated to form a 2D granular 
gas |T7H22]. Under mechanical agitation particles un¬ 
dergo a Brownian-like motion, but due to the dissipative 
nature of the collisions, the granular gas reaches a sta¬ 
tionary out-of-equilibrium state. Then immersed in a 
vertical external magnetic field B, the soft ferromagnetic 
spheres are magnetized and interact with each other as 
induced dipoles. When the magnetic field is increased, 
the granular gas solidifies into a phase composed of chains 
of a few particles in contact, similar to the labyrinthine 
phase observed with colloids m- In contrast to this 
colloid study which focuses on equilibrium states, the 
transition from gas to labyrinth is here clearly described, 
using an accurate particle tracking. Finally as a remark, 
we emphasize that the physical mechanisms at play in 
labyrinthine and stripe phases of interacting particles, 
differ from those in chain and cluster phases reported 
in some interacting granular gases despite vi¬ 

sual similarities. Indeed, these phases are composed of 
head-to-tail dipoles and were observed when attractive 
behavior is dominant at large scale, because permanent 
dipoles are considered [531 [3^ or because the hypothesis 
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FIG. 2. (color online) Top views of the system of particles. In (a) and (b), images of the spheres have been replaced by white 
disks of diameter a for better visualization, whereas (c) is from a direct image from the camera and (d) is the result of particle 
tracking, (a) Dissipative granular gas state at moderate B (80G). (b) Labyrinthine phase at high B (170G). The particles 
organize into an amorphous phase mostly composed of chains of a few particles. The region within the red square is enlarged in 
(c) and (d). (c) Thick circles are reflections of lighting on the spheres and appear smaller than the actual overlapping particle 
radii (thin circles), (d) Buckled chains (red or light gray for top particles, blue or dark gray for bottom particles) mostly 
compose the amorphous phase. 


of a quasi-2D system is not verified 

II. FROM A GRANULAR GAS TOWARDS A 
LABYRINTHINE SOLID PHASE 

The experimental setup is similar to the one used 
in [35]. Soft ferromagnetic spherical particles of diam¬ 
eter a = 1 mm are confined between two horizontal par¬ 
allel plates separated by a gap h = 1.5 a in order to form 
a monolayer as depicted in Fig.j^ Particles are vibrated 
vertically and are lightened by an annular LED array and 
imaged from the top by a fast camera through the trans¬ 
parent top plate. Particle center positions are tracked 
and their trajectories are reconstructed in the horizontal 
plane |3TJ . Interactions between particles are intro¬ 
duced by means of an external vertical magnetic field 
of amplitude B controlled by the experimentalist. Ad¬ 
ditional details on the experimental setup and protocol, 
and on the particle detection technique, are given in the 
Appendices. 

Magnetized soft ferromagnetic spheres behave as in¬ 
duced dipoles, whose magnetic moments are vertical and 
proportional to B. The interaction potential Um of a 
particle located at a distance r and a polar angle 6 from 
a second particle [35] reads in spherical coordinates |33| , 
with /To the vacuum permeability: 

^m{r,0) = — ^ (2cos^6'-sin^6') (1) 

16 fXQ ^ ' 

Two spheres in the same horizontal {9 = 7r/2) plane 
are thus repelling each other. Using this experimen¬ 
tal method, macroscopic transitions were observed in 
3D assemblies of magnetized soft-ferromagnetic parti¬ 
cles [301 Elj. Then, the number of particles per surface 
unit is expressed by a dimensionless parameter, the area 
fraction (j) = {Ntto^)/{4 S), with N the number of parti¬ 
cles tracked in the region of interest S. For a monolayer of 


particles, high enough magnetic field and moderate area 
fraction {(f> = 0.2), the 2D granular gas solidifies into 
a hexagonal crystal [35], whose melting has been found 
to follow the Kosterlitz-Thouless-Halperin-Nelson-Young 
(KTHNY) scenario [53]. 

Here (f) is increased to 0.5. For moderate values of the 
magnetic field B and continuous shaking, particles un¬ 
dergo a Brownian-like motion. At a given instant parti¬ 
cle positions are random [Fig. ^a)] and spheres exchange 
energy through dissipative collisions and magnetic inter¬ 
actions. We observe a 2D granular gas, whose properties 
are similar to those found at lower area fraction [26] . 
Then by increasing B further, we observe that despite 
the magnetic repulsion, small chains of two or three par¬ 
ticles in contact are starting to form in the bulk of agi¬ 
tated particles. We remark also that the motions of the 
particles belonging to these chains are considerably re¬ 
stricted compared to those of the free particles. At higher 
magnetic field, the quasi-totality of the particles are con¬ 
densed into these chains [Fig. [^b)]. At large-scale, the 
picture of the assembly of the system does not present 
an ordered structure. Thus, by increasing magnetic in¬ 
teractions, the system has been solidified in an amor¬ 
phous state. Labyrinthine patterns were indeed described 
as globally disordered stripe domains [33]. Due to the 
presence of chains, the particle assembly presents strik¬ 
ing similarities with the labyrinthine and stripe phases 
observed [T5] and numerically predicted [TUHH] for ex¬ 
ample for 2D colloidal systems under thermal agitation 
and with dipolar repulsive magnetic interactions. Here, 
a transition from a granular gas phase to a labyrinthine 
phase for a macroscopic, out-of-equilibrium system is ob¬ 
served. This transition can be also visualized by apply¬ 
ing a linearly increasing magnetic field from H = 0 G to 
B = 200 G (see Movie in Supplemental Material [34] ). 
First, pairs aggregate, then triplets, and so on homoge¬ 
neously across the cell, until nearly motionless chains of 
various lengths occupy the whole cell, isolating the few 
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remaining fluctuating particles from each other. We no¬ 
tice also that, starting from the labyrinthine phase, the 
inverse transition is observed when the magnetic held is 
decreased. This shows that the system does not present 
any noticeable hysteresis. 

By the means of an accurate particle tracking, chain 
morphology can be now quantitatively characterized. 
Chains are well separated due to the magnetic repulsion 
and can be thus considered as groups of more than one 
particle, according to a criterion of the contact distance. 
Moreover chains appear mainly as linear objects because, 
most of the time, a particle inside a chain is in contact 
with two neighboring particles. Nevertheless the relative 
orientations of the chains seem random. In the following, 
quantitative analysis of the small-scale structure will re¬ 
veal that chains correspond to a buckled state of particles 
in contact [Fig. [^c)]. A particle, once it is condensed in 
a chain, is either in contact with the top plate or with the 
bottom plate. Using slight differences of lighting for the 
two kinds of particles, our detection technique is able to 
provide the vertical position of the particles in the chains 
[as shown in virtual image Fig.j^d)], that is coded as up 
(red) or down (blue). It can be noticed that the parti¬ 
cles at the tips of the chains are for the most part down, 
revealing an effect of gravity. 

III. CHARACTERIZATION OF THE 
TRANSITION 

Let us now quantitatively characterize the transition 
from granular gas to the labyrinthine phase using the par¬ 
ticle tracking data. The magnetic potential energy Em 
and the horizontal kinetic energy Ec per particle can be 
now computed. From the interaction potential Um(j’-,&) 
defined in Eq. Q, the magnetic energy per particle is 
computed as the averaged summation over the pairs of 
the interaction potential of each pair: 

iVp iVp 

^rn = Um(j’ij,6ij) ( 2 ) 

^ i—1 j—i+1 

with Np the number of particles involved in the calcu¬ 
lation of Em, Uij and 9ij respectively the distance and 
the polar angle between the two particles i and j, and 
• the temporal average. The magnetic potential energy 
depends on the local configuration of the particles, and 
therefore it fluctuates in time. Its averaged value, Em, is 
found to be proportional to [Fig. [^a)]. The kinetic 
energy per particle is computed from velocity measure¬ 
ments: 

+ (vy^)) (3) 

where m denotes the particle mass, Vx (resp. Vy) the 
horizontal velocities in the x-direction (y-direction) and 
(•) an ensemble average. Ec is a measure of the agita¬ 
tion in the system. When the magnetic held B is in¬ 
creased, in a Hrst stage Ec grows [Fig. |^b)]. Repulsive 
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FIG. 3. (color online) Potential magnetic energy and kinetic 
energy per particle. Every marker corresponds to an indepen¬ 
dent experiment, (a) The potential magnetic energy is found 
to scale as (the dashed curve is a fit), (b) The kinetic 
energy, which measures agitation, plotted as a function of B. 
First, Ec increases due to the fluidizing effect of the magnetic 
interactions [26], then it suddenly drops towards zero at the 
onset of the solidification {B ~ 100 G). 


dipole-dipole interactions reduce the rate of dissipative 
collisions, which increases consequently Ec for a given 
shaking amplitude [21]. Once chains start to form, for 
B « 100 G, Ec drops signihcantly and nearly vanishes as 
the labyrinthine phase is formed for B « 150 G, illustrat¬ 
ing the solidihcation process. For higher magnetic ex¬ 
citation values, the labyrinthine phase becomes less and 
less mobile as interactions strengthen and fluctuations 
are restrained. 

Now, let us define the dimensionless control parame¬ 
ter e = Em/Ec [26|, which is depicted as a function of 
B in the inset of Fig. |^a). e provides a measure of the 
competition between distance interactions and kinetic ag¬ 
itation. By analogy with an order parameter, the frac¬ 
tion of particles condensed in the chains (p is computed 
as the ratio of the number of particles belonging to a 
group of more than one particle to the total number of 
particles tracked in the region of interest. By plotting 
(/? as a function of e, as shown in Fig. |^a), the transi¬ 
tion is well depicted. For e < 60, is nearly null in the 
granular gas phase, whereas for e > 1000, ip is slightly 
smaller than one, for the labyrinthine phase. The inter¬ 
mediate region of partial solidification corresponds to a 
coexistence zone between fluidized particles and parti¬ 
cles condensed in the chains. Let us emphasize that, al- 
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FIG. 4. (color online) Characterizing the transition. Every marker corresponds to an independent experiment, (a) The main 
plot shows the fraction of particles in chains (p, as a function of the control parameter e = Em/Ec. It varies from about zero 
in the dissipative state to one in the fully solidified state. The horizontal, black dashed line is placed as a guide. The inset 
shows e as a function of B. (b) Adimensionalized mean chain length A/a as a function of e, with A evaluated by computing the 
largest distance between particle centers inside a given chain. The value of A averaged over the chains increases continuously 
with e. If the averaging is weighted by the number of particles in the chains, the obtained values are larger by roughly one 
diameter unit, but the trend with e is similar, (c) Evolution of the radial pair correlation function g{r/a) for growing values 
of e (see arrow): e = 14.7 (red), e = 41.0 (purple), e = 68.9 (green), e = 127 (orange) and e = 649 (blue). As the transition 
occurs, the r = a peak of g{r/a) drops down to zero, showing that in-plane collisions probability vanishes for high values of e. 
In the meantime, another sharp peak grows at r/a = Xg/a = 0.91, which is the imprint of the buckled chains. 


though the simple criterion of the variation of p) captures 
well the transition from a granular gas to a labyrinthine 
phase, it should fail to distinguish a stripe phase from 
a labyrinthine phase. Topology and morphology would 
indeed have to be taken into account, like local orienta¬ 
tional properties. Several approaches were proposed to 
analyze or to model labyrinthine patterns, such as the 
introduction of a local wave vector [35], the computa¬ 
tion of the wrinkledness [36] . and the decomposition of 
the pattern into clusters of linear segments |33j . To our 
knowledge, the definition of an appropriate order param¬ 
eter for labyrinthine patterns remains an open question. 

Nonetheless, aiming at quantifying some of the 
directly observable morphological changes of the chains, 
we evaluate their mean length A as a function of e 
[Fig. gb)]. A is defined as the average over all chains 
of the largest distance between particle centers inside 
a given chain. Starting from 1.6 o at the formation 
of first chains, A seems to saturate for the highest 
values of e around 2.6 a. This suggests that compe¬ 
tition between growing chains could limit their extension. 

The pair correlation function, which is related to the 
probability of finding a particle center at a given distance 
from another particle center, provides information on the 
small-scale structure of the system. Before the transition, 
the particle assembly evolves from a purely dissipative 
gas to an effectively more elastic gas [3^. Therefore, as 
Fig. I^c) (red and purple lines) displays, the peak at the 
diameter value flattens whilst the effective elasticity rises 
up. From the onset of solidification (green, orange and 
blue lines), surprisingly an extremely sharp peak grows 
from zero at the distance value r = Xg ~ 0.91 a, which 
is smaller than the particle diameter. This peak, which 
would be impossible to observe in purely 2D systems of 


hard spheres, reveals the internal structure of the chains. 
Here the gap size is indeed large enough so that partial 
overlaps of particles are allowed, leading to the formation 
of buckled particle chains in which particles are in contact 
with the top or bottom plate [see schemat ic in Fig.[^ c)] . 
Geometrical calculations yielding Xg = \/2ha — h'^ , one 
finds h = 1.42 mm, which corresponds to the announced 
gap of 1.5 mm diminished by the the roughness of the 
bottom plate. Moreover in the labyrinthine phase (blue 
line), g{r) exhibits also a shorter peak at the position 
2Xg, from the aligned second neighbors, showing the 
presence of linear chains. Between these two peaks, a 
zero-probability at Xg ^ r < 1.2 a indicates the void 
spaces between the chains while the non-vanishing prob¬ 
ability for 1.2 a < r < 2xg stands for both non-aligned 
second neighbors in chains and particles from neighbor¬ 
ing chains. The 3D effects related to the gap size h are 
thus essential to describe the small-scale structure of the 
labyrinthine phase. Therefore, we discuss now how three- 
dimensionality can explain the stability of the chains at 
high enough area fraction. 


IV. CHAIN FORMATION MECHANISM 

At low area fraction (</) = 0.2 and lower), the stable 
state of the assembly of spheres in dipolar interaction was 
found to be a hexagonal crystal |26l [32] . Why does the 
hexagonal structure become now unstable at higher area 
fraction and how can we explain the formation of chains 
of particles in contact? In Fig. [^a), we plot the 3D mag¬ 
netic energetic landscape (in a vertical plane) for a sphere 
initially in the center of a hexagon of 6 neighbouring par¬ 
ticles, the projected horizontal distance between the par- 
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FIG. 5. (color online) (a) The potential energy landscape Em 
is computed for 0 = 0.5 and B — 200 G by averaging the pair 
potential Um [see Eq. Q] for six down neighbors forming 
a hexagon (four particles are out of the figure plane) and 
a central up particle moving along the x-axis (dashed line), 
(b) The proEle of Em is plotted along this trajectory for four 
values of (p (curves have been rescaled and shifted vertically 
for clarity). Circle mark depicts the initial central position 
and cross marks show the two contact positions. From these 
graphs, for cp = 0.2 and 0.35 hexagonal conhgurations are 
found stables, whereas for cp = 0.5 and 0.6, the central particle 
in presence of agitation should leave the position x = 0 to 
reach contact positions associated with the buckled chains. 


tides being geometrically given by d = a y 7r/(2-\/3(/)). 
Let us consider this central particle at x = 0 assumed up 
[dashed circle in Fig.j^a)] between six down neighbors all 
at a distance d, and investigate its potential energy when 
it moves from x = 0 to the contact position for several 
values of (p [Fig.j^b)]. Contact positions between spheres 
are local minima of potential energy, as the dipolar in¬ 
teraction aims to align spheres along a vertical axis. For 
(p = 0.2, the central position is an absolute minimum of 
energy, in agreement with the expected stability of the 
hexagonal lattice. In contrast when (p is increased, the en¬ 
ergy barrier decreases and Em{Q) augments relatively to 
ifm(±d). For (p = 0.5 contact positions become absolute 
minima, which can be reached by means of mechanical 
agitation. For (p = 0.6, the central position is not even a 
minimum anymore. Therefore for (p = 0.5 and above, we 
expect that hexagonal structure is unstable, leading to 
local structures of spheres in contact like chains, despite 
the isotropic dipolar repulsion in a purely 2D system. 
This qualitative model explains the small-scale attrac¬ 
tion leading to particle contacts, needed for the shaping 
and the stability of labyrinthine phases [5] . To improve 
the description, solid friction between the spheres and 


the top and bottom plates should be incorporated, as 
it may greatly stabilize the buckled chains. In thermal 
systems, similar predictions were obtained using Monte- 
Carlo simulations HOHIl]. In both these examples and in 
our system, the resultant of repulsive interactions of the 
assembly of particles over one acts as a magnetic pres¬ 
sure, favoring contact at high enough particle density. 
Additionally, we notice that buckled phases stabilized by 
pressure and friction can also appear in thin vibrated 
granular layers without magnetic interactions [37j if the 
density and gap size are sufficiently large [35]. Neverthe¬ 
less, in this last case, the structuring in separated chains 
is absent. 


V. CONCLUSION 


In this macroscopic and out-of-equilibrium model ex¬ 
periment, a labyrinthine phase is obtained by applying a 
magnetic field to a confined granular gas, by the means 
of externally controlled dipolar interactions. We describe 
and quantitatively characterize the transition from a gas- 
like phase towards a globally disordered solid phase. It 
appears as a three-dimensional effect occurring in a quasi- 
two-dimensional system. The parameters setting the con¬ 
finement, the gap h/a and the area fraction <p, are thus 
essential to explain the phase diagram of this granular 
medium, as it was also shown for colloidal and hard 
spheres monolayers [I2lll6l|39|. Although not presented 
here, after a fast increase of B, i.e. a magnetic quench, 
the labyrinthine phase exhibits a slow dynamics char¬ 
acterized by a slow evolution of its structural proper¬ 
ties [40]. This aging phenomenon should thus be com¬ 
pared to the slow dynamics of structural glasses |3Tti35] . 
In the case of labyrinthine domain patterns arising in 
continuous systems, the analogy with glasses was also re¬ 
ported in the analysis of their globally disordered struc¬ 
ture |33j and when studying their relaxation |36j . Fi¬ 
nally, whereas the structure of the phases obtained in 
this macroscopic experiment resembles those found at 
thermal equilibrium in Monte-Carlo simulations and col¬ 
loidal monolayers [Ulllllll], kinetics of the transition 
described here is intrinsically an out-of-equilibrium pro¬ 
cess, which deserves further studies. 
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APPENDIX A: POTENTIAL MAGNETIC 
ENERGY FOR TWO PARALLEL DIPOLES WITH 
FINITE MAGNETIC PERMEABILITY 


For two ferromagnetic spheres of identical diameter 
a immersed in an unidirectional vertical magnetic field 
of intensity B and separated by a distance r, the po¬ 
tential energy of magnetic interaction reads in spherical 
coordinates as the interaction of two vertical magnetic 

dipoles [55] : 


Umir,9) 


TT a® / fi — fj-o 

16 fio 


2 

(2 cos^ 


sin^ 0) 


(4) 

where 0 is the polar angle between the two dipoles, 
/io = 47rlO“^H.m“^ is the vacuum permeability and /r 
is the intrinsic magnetic permeability of the sphere ma¬ 
terial. The induced magnetic fields of the neighboring 
particles are assumed negligible in front of the external 
magnetic field. For free moving particles 9 is taken equal 
to 7r/2, whereas for particles belonging to chains 9 is com¬ 
puted from the measured vertical position of the parti¬ 
cles (top or bottom). The magnetic potential energy per 
particle is computed as an average of Um over inter¬ 
acting pairs of particles |26| and in the limit of large /r. 
This approximation holds for soft and linear ferromag¬ 
netic materials [45] . which is the case for our particles. 


APPENDIX B: GENERATION OF THE 
VIBRATED AND INTERAGTING SYSTEM OF 
PARTIGLES 

The particles are chromed alloy steel (AISI 52100) 
spheres of a diameter a = 1 mm and of a mass 
m = 4.07 X 10“^ g. Using a vibrating sample mag¬ 
netometer, the magnetization of one particle was mea¬ 
sured by V. Dupuis. The magnetic permeability /r 
verifies 122 < < oo, in the linear domain 

(—2000G < B < -1-2000 G) and the remnant magnetic 
field Br is below 12 G. The coercive field is small com¬ 
pared with the values of magnetic field B used in our 
experiments. Within this range of B, the response to 
magnetic excitation is linear. The square aluminum cell 
(side 9 cm long) containing these particles (see Fig. 
is vertically driven by an electromagnetic shaker. The 
forcing is sinusoidal at the frequency /o = 300 Hz and 
the root mean square (RMS) acceleration of vibration is 
fixed at F = 21.9 m.s“^ = 2.23 g for all experiments, with 
g the gravity acceleration. This value corresponds to the 
upper limit of the linear response domain of the granu¬ 
lar temperature Tg = Ec/m as a function of F [21] . The 
two coils are a Helmholtz pair current-controlled supplied 
generating a nearly homogeneous vertical magnetic field 


B across the cell (3% of variation measured). Immersed 
in this magnetic held, the particles are magnetized into 
induced dipoles vertically oriented (particle rotating ve¬ 
locity is negligible compared with the speed of the mag¬ 
netic domains rearrangements). 


APPENDIX G: PARTICLE DETECTION 

An annular light-emitting diode (LED) array above 
the cell produces a high-contrast circular signal on the 
chromed particles whose positions are recorded from 
above using a high-speed video camera at high resolu¬ 
tion (1152 X 1152 pixels at 780 Hz). The region of interest 
S is 5.7cmx5.7cm around the cell center (see Fig. [^. 
The particle diameter is about 20 pixels. For individ¬ 
ual particle detection, we used a convolution-based least- 
square htting routine [211146] completed by an intensity- 
weighted center detection algorithm (accuracy estimated 
less than 0.3pixel). Particle trajectories were recon¬ 
structed using a tracking algorithm EZlllT]. 


APPENDIX D: EXPERIMENTAL PROTOGOL 

It is fully automated for the sake of robustness. Ev¬ 
ery single experiment is non-correlated with the others. 
The amplihers of the electromagnetic shaker and of the 
Helmholtz coils are computer-controlled via a data acqui¬ 
sition card. The experimental protocol routine is written 
in Matlab. It also proceeds to the dialogue with the cam¬ 
era, i.e., configuring and starting the video recordings, as 
well as to the recording of the data from the accelerom¬ 
eter and the Hall effect sensor. All experiments are per¬ 
formed according to the following protocol. Eirst the 
shaking is activated (F = 21.9 m.s“^) while the magnetic 
field remains zero. The magnetic field is then linearly 
increased (rising rate a, = dB/dt is kept hxed for all ex¬ 
periments) up to its higher plateau value B. Afterwards, 
a waiting time is respected prior to proceeding to the 
recordings. It is chosen along with the recording time 
length so as to reach the chosen mean aging time Tw¬ 
in all the experiments presented here, ag = lG.s“^, 
Tw = 30 s and recordings last at least 2 s. Note that 
these two parameters Uq and r^, have a noticeable influ¬ 
ence on the nature of the labyrinthine state reached for 
high values of H, implying that a slow dynamics is at 
play. 

For the 5000 particles introduced in the experiments, 
the area fraction evaluated on the cell is equal to 0.485. 
However, within the region of interest S', as the bound¬ 
aries are not repulsive [55], (j) decreases from 0.58 to 0.46 
with B until H « 80 G. From the appearance of the first 
chains, (j) remains nearly constant. 
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